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ABSTRACT 


An analysis of steady state operation of a vertically 
oriented, variable conductance heat pipe is presented. The 
effects of binary mass diffusion, axial pipe wall conduction, 
and gravitational effects, caused by the difference in 
molecular weight of the non-condensible gas and the working 
fluid are incorporated. Analytical expressions for the con- 
servation of mass, momentum, and energy are combined along 
with equations of state to describe steady state operation 
of the heat pipe. These expressions are combined to form 
a system of three differential equations with three unknowns; 
working fluid vapor velocity, mass fraction, and temperature. 
These equations are nondimensionalized and an iteration scheme 


for numerical solution of the equations is presented. 
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axial conduction heat flux 
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working fluid vapor saturation density as a function 
of temperature 


absolute temperature 
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A. BACKGROUND 

A heat pipe is a closed, evacuated chamber that has the 
inside walls lined with a wick structure, saturated with a 
volatile working fluid. The steady state operation involves 
vapor heat transfer and capillary action. A comprehensive 
coverage of the theory and operation of heat pipes can be 
found in Dunn and Reay [Ref. 1]. Briefly, the operation of 
a vertically oriented heat pipe will be discussed. The heat 
pipe has three sections: an evaporator section, an adiabatic 
section, and a condenser section. The heat pipe with the 
three sections marked is shown in Figure 1.1. The lower evapora- 
tor end of the heat pipe is heated, causing the working liquid 
in the wick to vaporize. This phase change causes the working 
fluid to absorb latent heat and the vapor pressure increases 
in the evaporator end. The working fluid is drawn away from 
the wick, causing the liquid pressure to decrease. The vapor 
moves to the cooler condenser where it is condensed and the 
working fluid gives up its latent heat. The condensate is 
deposited on the wick which causes the liquid level to in- 
crease thus increasing the liquid pressure. Capillary forces 
take the higher pressure liquid from the condenser to the lower 
pressure liquid end of the evaporator. The capillary forces 


are assisted by gravitational forces when the heat pipe is 
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Poaemeaitly OLremred aS Shown im Figure I.1. It is important to 
note that for heat pipes the evaporator position is not re- 
stricted, but since the experimental data and numerical model 
use a vertically oriented heat pipe, this orientation will be 
used throughout this thesis. 

The great advantage of heat pipes is that the vaporization 
and condensation process allow for high heat transfer rates 
with small temperature differences. Two other advantages are 
that the basic heat pipe has no moving parts and can be used 
in many geometrical configurations. The heat pipe has limi- 
tations and again, a comprehensive presentation of these limi- 
tations can be found in Dunn and Reay [Ref. 1]. Briefly, the 


five heat pipe limitations are as follows: 


(a) The viscous limit is basically the low temperature limit. 


This would be the minimum temperature for the heat 
pipe to begin operation. Generally, the viscous limit 
is only important during start up of some heat pipes, 
1.e., liguid metal pipes. 


(b) The sonic limit describes the limit of the vapor flow 
velocity at the heat pipe evaporator exit. The re- 
sulting choked flow condition limits the axial heat 
flux capacity of the heat pipe. 


(c) The entrainment limitation is due to the removing of 
the liquid before it can return to the evaporator by 
vapor flowing at high velocity. 


fay The wicking or capillary limit reflects the fact that 
each combination of wick and fluid has a maximum 
capillary head. As stated earlier, the capillary head 
must increase with the liquid and vapor pressure drops 
and can do so only up to the maximum capillary head. 


feos ihe last limit is the boiling limit. If temperatures 
in the evaporator get too high the working fluid will 
begin to boil within the wick. This will disrupt liquid 
flow to the evaporator, causing the heat pipe to stop 
working. 
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These five heat pipe limitations vary in importance, depending 
on the temperature and corresponding axial heat flux operating 
point. Since at different temperatures one can encounter 
different heat pipe controlling limitations, @t seeaimpomecnme 


to have a basic understanding of all five imeracren—- 


B. VARIABLE CONDUCTANCE HEAT PIPE 
I. Basic Operatron 

Variable conductance heat pipes are described in detail 
in Marcus [Ref. 2]. As described in the background sectii@mpurs 
change in the heat load will automatically change the operating 
temperature and pressure of the heat pipe, given a constant sink 
temperature. This is an undesirable quality for many situations 
where it is desired to have a stable source temperature under 
conditions of varying heat load or sink temperatures. An 
example of this case would be the cooling of sensitive elec- 
tronic equipment. Basically, a noncondensible gas is introduced 
into the heat pipe along with the working fluid, making rhe 
heat pipe become a passively controlled device. The nonconden- 
sible gas occupies a certain portion of the condenser secticny 
depending on operating conditions. Increasing operating tenpewes 
ture in turn increases the vapor pressure of the working fluid 
which compresses the noncondensible gas into a smaller volume. 
This provides a greater active condenser area, thus allowing 
a more stable source temperature than a conventional heat pipe. 


When the operating temperature decreases, working fluid vapor 


IZ 


Pueactme talls and Che nonmcondensible gas expands, blocking a 
larger portion of the condenser. The result is a passively 
controlled variable condenser area, which in turn makes a 
more stable source temperature. 
ae eee wecOnguctance Heat eLpe mheory 

Marcus [Ref. 2] describes a flat front theory anda 
feet use £ront theory” for variable conductance heat pipes. 
Flat front theory assumes the noncondensible gas is in the 
form of a plug with sharply defined vapor-gas interface which 
is perpendicular to the axis of the pipe. Since this model 
neglected the axial conduction in the heat pipe wall and diffu- 
Sion of vapor into the stagnant gas plug, experimental results 
(Ref. 2] have shown this model to be inaccurate. The early 
model of flat front theory provided initial insight into the 
Operating characteristics of gas loaded heat pipes. With the 
incorporation of binary diffusion between the working nea BIBL 
Vapor and the noncondensible gas and axial conduction in the 
pipe wall, Marcus' [Ref. 2] ‘diffuse front theory' establishes 
a model of gas loaded heat pipes which provides better simi- 
larity between computer program and laboratory results. 

Gravitational effects on the vapor-gas interface are not 
fmermuded in flat front theory or the ‘diffuse front theory. ' 
Gravitational effects can greatly distort the axial temperature 
profile when the working fluid and noncondensible gas have 
Significantly different molecular weights. This fact 1s demon- 


strated by experimental results reported by Kelleher [Ref. 3]. 


3 


Batts [Ref. 4] conducted variable conductance heat pipe experi- 
ments using working fluid/noncondensible gas combinations of 
methanol/helium and methanol/krypton. Kane [Ref. 5] developed 
an analytical model of variable conductance heat pipes which 
includes the effects of binary mass diffusion, axial pipe-wall 
conduction, and gravitational effects which exist at the vapor- 


gas interface region as shown by Batts [Ref. 4] experiments. 
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If. FORMULATION OF AN ANALYTICAL MODEL 


A. PROBLEM STATEMENT 

The purpose of this chapter is to formulate the equations 
describing the steady state operation of a vertically oriented, 
@as loaded, variable conductance heat pipe. As in Figure l1.l, 
the condenser is above the evaporator for this orientation. 
This analysis includes the effects of binary mass diffusion, 
axial pipe wall conduction, and gravitational effects due to 
the difference in molecular weight between the noncondensible 


gas and the working fluid. 


Eee fnooUMPTIONS 
The following assumptions are used to develop the system 
of equations: 


l. Heat, mass and momentum transport processes are treated 
as One-dimensional in the axial direction within the 
vapor space of a vertically oriented heat pipe. 
Kelleher [Ref. 3] shows that in any other orientation, 
gravitational effects cause highly three-dimensional 
Spatial variations when the molecular weights of the 
working fluid and noncondensible gas differ substantially. 


2. The noncondensible gas acts as a perfect gas. The 
mixture of noncondensible gas and working fluid are 
EEeatead as thermodynamically independent substances. 


3. Because of direct contact between the vapor and the 
liquid saturated wick, the working fluid vapor is at 
saturation conditions at all times. 


4. The thermal resistance of the wick-wall structure, in 
the radial direction, is negligible in comparison with 
the thermal resistance due to convection outside the 
pipe. 


aes. 


Co. ANALYS US 

Kane [Ref. 5] presents a proposed solution for this prob- 
lem. The following is a significant change to his proposed 
solution with modifications being noted as they occur in 
this section. Variables with the subscript ‘a' apply tome 
working fluid and variables with the subscript 'b' apply to 
the noncondensible gas. 

Utilizing the assumption that the noncondensible gas is 
a perfect gas, the equation for noncondensible gas pressure 
1S: 


P = ak 


b PL i, al (Ii 


The principles of mixtures of thermodynamically independent 


substances, which are found in Reynolds and Perkins [Ref. 6], 


are used to define the total density and the total pressure. 


Cr -— e) aa oe (rT ae 
P = P +P (TT 


The definition of mass fraction is used to form the equation 


for the working fluid vapor mass fraction. 


W) = — (TI. 4) 
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The solution, utilizing the conservation of mass, is 
based on a differential element representing a section of 
the heat pipe condenser, Figure 2.1. The conservation of 


mass of species ‘a' for the differential element yields: 


_ d : 
O VA = On We + qx (Pag) AXA, ss m_P dx 


This equation sets the mass of species ‘a' into the control 
volume equal to the mass of species ‘a' out, including a 

term to account for the removal of species ‘a' by condensation. 
The process of condensation is accounted for by multiplying 
the condensation mass flux times the area, defined by the 
vapor space perimeter times the length, dx, of the differen- 
tial element. After basic algebra, the equation for conser- 


vation of mass of species ‘'a' becomes: 


Peeesing Fick's law of binary diffusion [Ref.7], and solving 


imep@ PAY! the resulting formula is: 


The final form for the conservation of mass for species '‘a' 


becomes: 


= 0. (Die) 


SSeW] JO UOTFEAATOSUOD TOF JuowSeTY TetjuessFjtqd “T°? Se ee 


1 24 a 
Bp (ea s+ TA Tg Bape ae eae d bps 


6 On Eos 
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Since species 'b' is a noncondensible gas, the condensation 
mass flux of species 'b' equals zero. The resulting equation 


for the conservation of mass for species 'b' is: 
S—(w,0V) - S-(pp =) = 0 
Goan. asx ax 


The mass continuity of the mixture is derived by adding 
| 


the conservation of mass of species ‘'a' with the conservation 


of mass of species '‘'b'. 


a d d GW. uy, 
eee ab!) da Ody * “ax? > © 
Because wW_+W, 1S equal to one, and Sern +wW,) = a) = 0 
a b ‘ axa b ax . 


the final result of the mass continuity of the mixture is: 


: Ee d 
m A a qx fev) = QO. (i146) 


a 

Figure 2.2 is a differential element of the vapor space in 
the condenser and is used to derive the conservation of linear 
momentum equation for the mixture of noncondensible gas and 
working fluid vapor. The conservation of linear momentum for 
the mixture of noncondensible gas and working fluid vapor is 
derived by setting the summation of external forces acting on 
the control volume in the x direction equal to the net efflux 


of linear momentum. 


Mm 
ty 
l 


> > > 
7 f V(pV dA) 
GS7 x 


aie, 
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The summation of external forces acting in the x direction are: 
YF = Ad d dp A 
= 9 9A, dx tL a x p A, 


These terms represent the force due to gravity acting on the 
control volume minus the forces due to wall shear stress and 
pressure acting on the control volume. The wall shear stress 


eam be calculated by: 


ay, 8 4 

le Re As V 
where X is a dimensionless friction factor. If developed 
laminar flow is assumed the Darcy friction factor, i = =<, 


could be used with the Reynolds number defined as 


pe VD, F 


Re = 
u 


where yu is the viscosity of methanol vapor, V is the velocity 
Of the vapor and dD, is the diameter of the vapor space. The 
moms to represent the net efflux of linear momentum are: 


° Z 
VimiP ax + (ov +d(pv) (vtdv) ) A, — np A. 


> > > 
f V (ovdA) 
Gs ,X 


Vote wet wae ovay + A vd(ov) 
aav Ce: C 


Zell 


By setting the summation of external forces in the x directa 
equal to the net efflux of linear momentum and combining like 
terms, the result is: 
Ie m 1 
— =7eg te Sov") = — = aa (II. 7) 
6 a Cc V 
The above equation represents conservation of linear momentum 
for the mixture of noncondensible gas and working fluid vapor. 
Kane [Ref. 5] assumes zero velocity for the noncondensible gas 
and derives an equation for conservation of linear momentum 
of species ‘'b.' The species 'b' conservation of linear momen- 
tum equation is not necessary to derive the model in this 
analysis. Since the conservation of linear momentum for the 
mixture of the noncondensible gas and working fluid vapor is 
the more consistent approach, Equation (11.7) is used througme 
out this thesis. 

Figure 2.3 is a differential element of the condenser wall 
and represents the wick and pipe wall combination as one ele- 
ment. The wick-wall combination is treated like a fin and is 
used to derive the equation for the conservation of energy. 
Heat transfer is represented by the latent heat deposited on 
the 'fin' by condensation of the working fluid vapor plus the 
heat transfer from Fourier's law of conduction enteringgaa. 


15a @ are 


in a fg 


De 


fg 


P ax 
Vv 


Ptoqure 2,3. 








C pdx 
X 
hP_(T = clos 
dx 
S See Oat 
=< Ha. Gue "GFE. aaa oo 


Differential Element for Energy Balance 
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Heat transfer out 1s represented by Fourier's law of conduction 
leaving the 'fin' plus the energy which is convected away fnem 


Pheer 2 


Q =~ 0. +:dQ_ + h P. (lec 


The equation for conservation of energy is derived by setting 


heat transfer in equal to heat transfer OWE of Ene hae 


2 hP h 
adat S co Vie = 

2 > eae tele) A ee Ti 
Gix p p 


The experimental data from Batts [Ref. 4] is used to 
correlate this analysis. The working fluid is methanol, 
CH,0H. The working fluid vapor in the heat pipe is in direct 
contact with the working fluid liquid an the wickeand sie 
therefore considered to be a saturated vapor at the local 
temperature of the wick-wall combination. Eubank [Ref. 8] 
tabulates data for the saturation properties of methanol. 
These tables are used to generate functions of local temperature 
for partial pressure, partial density and heat of vaporization 
of saturated methanol vapor. 


(a) The partial pressure of saturated methanol vapor is a 


function, Fa, of local temperature. 


P =) ae) em: 
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ioe parelal density of saturated working fluid vapor 


is a function, Sa, of local temperature. 


ae = Sa(T) Clie 
te)” The working fluid heat Of VapoOrization is a function, 
Bre! of local temperature. 
Dee = Beg OE! ( Lae) 


In summary, the eleven equations which are stated in this 
chapter form the foundation for the analysis of a vertically 
oriented, variable conductance heat pipe. The remaining 
chapters derive a numerical model and iteration scheme based 


on the eleven equations from this chapter. 
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Ill. NUMERICAL MODEL FORMULATION 


A. VARIABLE ELIMINATION 
The system of eleven equations can ultimately be reduced 
to a system of three equations by eliminating variables. 
Equation (II.9) is eliminated by using the equation of state, 
Pa Cr).2% “kor Pi The partial density of the saturated working 
flue vapor, - (Equation II.10) is eliminated by using Sage 
ee Equation (II.11), is used to eliminate Des The non- 
condensible gas pressure, Pie is eliminated by using the perfect 
gas equation of state (Equation II.1). The treatment of womeeng 
fluid vapor and noncondensible gas as thermodynamically inde- 
pendent substances gives total pressure as the sum of partial 
pressures and total density as the sum of partial densities. 
Total pressure is eliminated by using the sum of \partial pres- 
sures, Equation (II.3). The noncondensible gas density is 
eliminated by solving Equatwton (il. 2) eee On and replacing 01, 
Well > —O.- The definition for mass fraction of the working 
fluid vapor, Equation (II.4), is used to replace total densijae 
After making the preceding seven substitutions, four equa- 

tions with four unknowns are left. The four remaining depen 
dent variables which are functions of axial position varc- 

(a) mo working fluid condensation flux alan eee sec] 

(b)  V, velocity of the mixture {ft/sec 

(ec) Wor working fluid mass fraction [0 < oes 1] 


(d) TT) temperatiaes | Rie 
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The four remaining equations are: 








eS dad ,Sa 
(a) mM. a a AS a V) = Q il ae 
Cc a 
P dw 
V d de Sab a _ 
(b) Jes AS a5 qx (Sa V) - rea a wate = 0 i id lg 2) 
d " 
a b 
P m p 
Be Sa Sy Scie 2 ly eee | Wao 
= Sa g) Ae ) ea De Gert . 3) 
a a Cc a Cun 
2 P Piers 
dar S _ fgv 2 
(d) 7 ran law Gl SURE es a mo = Q (III.4) 
obs Pp p 


Furthermore, Equation (III.1) can be eliminated by solving this 
equation for working fluid condensation flux, ma, and substi- 
tuting that into the three remaining equations. 

The resulting three equations with dependent variables of 


velocity, working fluid mass fraction, and temperature are: 





dw 
di {Sa d d ,Sa a ss, 
a a 
dFa d aL R 
(b) rane 1 qx (Sa (— -1) M Te) 
eee = Sees) 4 Fa( 82 y) Pvt WP Gee 6) 
7 Ww dx ‘WwW W ax ‘WwW ADL : 
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H pen 
(ec) —3 ~ 9, HCE St) - =35 ee = (Ill oe 
Pp a 


This set of three equations and three unknowns describe the 
Operation of a vertically oriented, gas loaded, variable conduc- 


tance heat pipe. 


B. BOUNDARY CONDITIONS 

These equations describe the thermodynamics of the heat 
pipe condenser. The independent variable, the distance along 
the heat pipe, is measured with the origin at the condenser 
end farthest from the evaporator. The total condenser length 
is defined as L. Based on the system of three equations, two 
boundary conditions are required for temperature and mass frac- 
tion and one boundary condition is required for velocmae 

From experimental data, the temperature at the far end of 
the condenser and in the evaporator are known. The boundary 


conditions for the dependent variable of temperacumemamer 


when x = 0, T = T_, where T_ equals the temperature at the 
Cc c 
far end of the condenser: 
when x = L, T = T_, where T_ equals the temperature at 


the end of the condenser closest to 
the evaporator. 

The highest velocity of the working fluid vapor is reached 
in the evaporator exit. This velocity can be calculatedytuem 
the known heat input to the heat pipe. The vapor velocity in 


the evaporator, Veo is equal to the heat transfer rate into 
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the evaporator divided by the product of local density of the 
working fluid vapor, cross sectional area of the vapor space, 
and the local heat of vaporization of the working fluid. Due 
to the physical boundary at the end of the pipe, where x = 0, 
the velocity at the farthest end of the condenser away from the 
evaporator is equal to zero. The boundary conditions for the 


dependent variable of velocity are: 


lI 
o) 
S 
oy 
a) 
Qu 


when x = 0, V 


when x = L, V= Ves = ae a 
a Ge tg 
e e 


Based on the system of three equations, only one of the above 
boundary conditions is required. The boundary condition at the 
Origin is used for the iteration scheme and ve 1s used to 
normalize velocity. 

The temperature and pressure at the farthest end of the 
condenser from the evaporator are known from laboratory data. 
Assuming thermodynamically independent substances, the total 
pressure is the sum of partial pressures and the total density 
is the sum of partial densities. The saturation pressure of 
the working fluid vapor, Poe Could be determined from tabulated 
data by entering the saturation pressure table with condenser 
temperature. An example of this table can be found in Eubank 
(Ref. 8] for methanol saturated vapor. The noncondensible gas 


pressure, P,, can be calculated by subtracting the working 


Ze 


fluid vapor pressure from the total pressure. With this dace 
the initial mass fractions can be established by using the 


following procedure from Reynolds and Perkins [Ref. 6]: 


P 
i = 
3 eOtal 
x es “b 
b ee eal 
x_M = ee 2 for species '‘a' 
aa Lbmoie mixture 
= Lbm ae 
Xbb Lbmole mixture FOL «SRG ye oaaag 
es 
U) —_ es 
ag nt eGR P 
Xp 
GU) = —— a 
Po Lees et 


At the end of the condenser closest to the evaporator, it is 
assumed that the concentration of noncondensible gas is zero 
SO; that at enusspoint Wo is equal to one. The boundary condi 
tions for the dependent variable of mass fraction for the 


working fluid vapor are: 


when x 


lI 
© 
= 

lI 
= 
a) 
= 
Qu 


when x 


lI 
be 
= 

lI 
— 
© 
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C. EQUATION NONDIMENSTONALIZATION 
TO generalize the equations and facilitate their solution 
it iS convenient to use dimensionless variables defined as 
follows. 
(a) The independent variable, condenser length, is defined 


Dy: 


= 
I 
yx 


where x is the distance along the heat pipe, measured from 
the condenser end farthest from the evaporator, and L is the 
total condenser length. 


(b) A dimensionless temperature can be defined by: 


where T is the local temperature, oe is the evaporator tempera- 
mie and T is the ambient temperature. In this analysis, the 
first and second derivative of temperature with respect to dis- 


tance are required. 


aT _ aT d0 dn _ “e “ede 
Sipe Ge an ax L dn 
and 
a*p _ a at, _ a te» ao,dn 
Bee dx “dx dn L cid x 
= ave d¢ 
te bac 


ol 


(c) The second dependent variable, working flgge@mass 
fraction, is dimensionless by definition and its numerical 
value is greater than or equal to zero and less than or equal 
to one. The first derivative of working fluid mass fractiuem 


with respect to distance is: 


dw. : dw. an _ 1 dw 
tian OF ccna L dn 
(d) Velocity of the working fluid vapor is made dimension- 


less by normalizing local vapor velocity with the evaporator 
Vapor velocity, Ve° The nondimensional velocity and the first 


derivative of velocity with respect to distance are then 











defined by: 
V 
ee 
V V - 
e 
V 
dx av scm. dix esta 
(e) In order to complete nondimensionalizing process, jeme 


following other dimensionless variables are used. 


(i) The working fluid vapor density is normalized with 


respect to the evaporator vapor density. 


SZ 


(ii) The working fluid vapor pressure is compared 


to the evaporator vapor pressure. 


Celene working Llugd@heat ef VapOrization 1S normalized 
with respect to the condenser working fluid vapor heat of 


condensation. 


fg 
k = 
eg H 


(iv) The working fluid vapor viscosity is normalized 


with respect to the evaporator vapor viscosity. 


After substitution of nondimensional terms, the three 
equations which describe operation of a vertical variable 


conductance heat pipe become: 


ad ,Sa*-Vv* d k : 
eo - ana, tap (Sa*Vv*) - —__(—_ ——_ —_— =O (Iles) 


The dimensionless coefficient from Equation (III.8) is 


D 


Nl =: 
= 


So 











Fae dFa* SO d pee i 
(b) — 2 oh wh “er an oo romanian! (9 ti Te 
sav, i a e 
bgusa® d Sa*y*? » GC Pasa 
> wl dan uw. ° \ aye 
Ve a a a 
Pp 8H L 
Tepe V7 Ls ee (TL i aay 


The dimensionless coefficients from Equation (III.9) are 











Fa, R(T_-T,) ae 
N2 = Ao N3 = =, cee N4 = ip ae 
Sa V May e a 
ee be 
Pee <S eee 
oe) Se) = V & 
N5 = y and N6 AbD Sav 
V GC Vv e e 
e 
Z Z 
a P AL H A Sa V_L 
(c) : et ees 9 Ene Cn ne rye 
ae K A, K A, (T, -T,) fg 
*KUT* 
x = (22 = eo (TTT an 
ig Ww. 


The two nondimensional terms of interest from Equation (III.10) 


are 


N7 = 


Numerical evaluation of the dimensionless CoOectificieneo ee 


required to determine the relative significance Of Gach tem 
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Two experimental runs from Batts [Ref. 4] are used to evalu- 

ate the dimensionless coefficients. A low powered run 

(Q = 22.9 watts) with methanol as the working fluid, and krypton 
as the noncondensible gas, is chosen. This run demonstrates 

an example of the heavier molecular weight noncondensible gas 
diffusing through a large portion of the condenser. A second 
higher powered run (Q = 70 watts) with methanol as the working 
fluid, and helium as the noncondensible gas, is chosen because 
it demonstrates an example of the lighter molecular weight 
noncondensible gas occupying a small portion of the condenser. 


The convection heat transfer coefficient for the heat pipe is 


Bew 


set egual to 1.0 oa aaa representing free convection. A 
hte. ° R 

system thermal conductivity of approximately 4.56 =e 

ier Gee ak 


results from a parallel circuit combination of conductive 
thermal resistances for the stainless steel pipe wall and 
methanol filled stainless steel wire mesh wick. Kane [Ref. 5] 
calculated representative binary mass diffusion coefficients 
of 0.4 ft-/hr for krypton-methanol and 1.9 ft*/hr Pore ne Li um= 
methanol. These values for binary mass diffusion are also 
used in this analysis. The experimental values that are used 
to calculate the dimensionless constants are listed in Table 
I. The values of the dimensionless constants for the two 
experimental runs are listed in Table II. 

Based on the tabulated numerical results, Equation (III.9) 


is modified to: 
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TABLE I 


NUMERICAL DATA 


Constants fOr bern cunsc : 


R = 1545.33 teen, Lomole °R Gas constant 

Guo a2 ft jae Acceleration of gravity 
Ln = 3S 5 ofc Length of condenser 

A, =~ OLSS fer Cross sectional area of 


vapor space 


) ere 01 0 em a een Axial cross sectional area 
P of wick-wall combination 

D. =" 15917 get Diameter of vapor space 

Pe =. 524) Et Periphery of outside of 
pipe 

Bae POO Oe tie Periphery of vapor space 

bh =. lO Bev eo or Convection heat transfer 
coefficient 

Ke] 4455965) beu/7 pre ork or Thermal conductivii, 


Low powered krypton/methanol 


QO = 22.9 wates = 78.1343 Btuyhe 


To =)59 22 
e 

T= 54372 
& 

T= 53608 


ey ee fine 
M = 83.8 Lbm/Lbmole 
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TABLE, f (CONTINUED) 


we =O al 3395 fe, SEC 
fis, = IU a ee 7 ee 
Heg,c =o. 2 Be / Lbm 


ep 
5) 
ll 


0.054871 Lbm/ft> 


g 


B= DOA See oLOm* bt Sec/ate 


High powered helium/methanol 


Q = 70 watts = 238.84 Btu/hr 


T = 640°R 
e 

T = 540°R 
Cc 

ie = 336°R 


Do= 1.9 ft*/hr 


My, = 4.003 Lbm/Lbmole 
Wi. = .050485 ft/sec 
ectg, 
Fa, =~eo. 2 bot /in 
Saat = 501.11 Btu/Lbm 
Ba = 0.14446 Lbm/ft~ 


2 


u, = 243.993 x10” Lbf Beer te 


a7 


TABLE. 


DIMENS TONLESS SCONS TAs 








Low Powered Krypton/Methanol High Powered Helium/Methanol 
a | ae -4 >” = -3 
Nie ihe 1 67 2 oo ME = a a 2.9869 x10 
e e 
Fa Fa 
N2 = f= 506.355 toe 10S f= 359.166 x 10° 
Sa V Sa_V 
e ee 
Bias) 6 Ra Fe 
Ne. = 5 =o192458)1.x 10 N3 = ———~ = £507.224 x10 
M, V M v7 
b ‘e be 
ine TO 
Na — 7, oT = yale N47 _J= 7 oT = 5 eee 
S 00 e CO 
NG. ao =o ~=eresoeaa N5 = ES) 2 Saeiieme 
2 2 
V V 
e e 
co Bu. E ee el, Wb 
N6 = rN D Sa V = 1A2086 N6 = AD. Sau = 4.7642 
cv ee V e 
Pp pare Pp AL’ 
ee z = eS 4 Sis N7 = a aaa 134.3552 
p p 
a y NS Y= en ee H A Sa VL 
Ng = EOG4:Cascc e e Ng = Gy, © ae => & 
(oy ie : K. Ae =e 
| oe = as 
= "53,302 = / LOOe as 
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Srdra~ IL 00 
——— + ———~——- [[—— (Sa* (— -1) (9 +=———_) ] 
fame sO" Move an Ne Be tes 
ee be 
* 
= cae reales) 
pe ci 
V a 
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The viscous term and the two inertia terms are neglected be- 
cause their order of magnitude is LO compared to TO and Lo 
for the remaining terms. The nondimensional equations (III.8), 


(I1I.10), and (III.11) are used to develop the analytical model 


Somene vertically oriented, variable conductance heat pipe. 


D. FORMULATION OF THE NUMERICAL MODEL 
Equation (III.8) is integrated and is solved for the first 
derivative of the working fluid mass fraction. Integrating 


both sides of Equation (III.8) gives: 








dw 
KV7* x 
ee oa Sl; Ciena) 
WW) dn 
a a 
The lower endpoint is used (n = 0) to evaluate the constant 


of integration. When n = O, V* at the end of the pipe equals 


Zemo. The constant of integration then becomes: 
0 dw 
o> N= Gunn ey 
Sa dn 
e 0 


By replacing the value for the constant of integration and 


Semving for dw /dn, EheseideteeverstOn Of ebquation (I1i.12)° 1s; 


So 


x 
2 a — a “a + —s —o a (IIT ae 
Nl Sa* Sa. dn 0 ° 14) 


The derivative in Equation (III.11}) is expanded to gets 

















following. 
6 pee 
an (88 Goo (6 +N4)) 
a 
dw 
2) edo ek pee a 
= an = Ip (ce +N4) + Sa al rare +N4) 
a W) 
a 
cae (=) (TIT ee 
W) dn 
a 
* 
The chain rule is applied to = 
ole 
dSa* _ dSa* dSa ay do = e »)dsa d6- (IITI.16) 
dn dSa dT dé dn Sa ela aoliq : 


* 
Applying the same procedure to oe , the result is: 





f= 
* CO 
SE = ) dra dé (IIT. lee 





By combining terms and solving for a the final form ¢en 


Squation: “Clie eed s. 
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(eumemelne COoectficient of a is set equal to G(wi,6). 


eT Se! dFa N3(T-T,) dSa 


il ; 
a = Fa ee Sa afc ot) he NS) 
e eG a 


rn N3 (Sa*) (= -1) (III.18) 
a 


(b) SOLVING fOr, a the result 1S 


dw 
7 Neipae N3Sa* a 
aT = a aS 2 “ae +N4)]/G(w_ 8) (elie os) 


a 





Equation ({I11.10) is integrated and the result is: 


Nn KUTk 
do f 6dn - N8 Hk sa*v* = 


CATT 220) 
0 tg Ya 


Meang the boundary condition of n O where V* = 0.0, the 


Semstant Of integration is: 


= ° (Tie i) 


dme final version of Equation (III.10) which is solved for V* 


1S: 


dé dé 
V fan Wo N7 (f6dn)w, ean Pedi > ae ae. (TEba272) 


4] 


The final system of three equations with nondimensional 


Varrab les as 














dw, dw. 
(a) Equateiren EE ey es as 
* 
du. _ vt V*w : W, Pp dw, 
a x Cy “AR 
dn N1 NL So Sa. dn 0 
dé oo 
(bo sEQUaE TOM ai ieisses jee ay ao ene 








dw 
clea NoSatee Nl sSare a 
ae [ > a. 5 an (8+N4) 1/6 (w, 8) 
a W) 
a 
where 
N2 - - 
Gieo = ten te) ee ee Te! SEE ea) (aad) 
a Fa el Sa ik Wo 
1 
+ N3(Sa*) (—— -1) 
a 
(c) Equation (III.22), v*(w ae Se ch 
q e , a’ "+ an’ an 
dé dg 
* = eee = = beudiody * * 
V ain Ws N7 (f8dn)w, W. an )/ (NSHE GS ) 


The boundary conditions for this set of equations abewed 


POW LOws Me 
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(a) For nondimensional temperature: 


ta 
Nn = 0.0 9 => 7 = oe 
e co 
fe) =e oo = 810 
dé; _ 
dn 0 


The constant a is not known but, based on laboratory data, can 
be estimated and through an iteration process the correct value 
can be obtained. A proposed iteration scheme is presented in 


Chapter IV. 


(b) For working fluid vapor mass fraction: 
ae — 0.0 a es “Ste 
n = 1.0 ee 1.0 
dw. 
ees = 8 


The constant § is not known, but can also be estimated for an 
gnitial value and, through an iteration scheme, a correct value 


can be obtained. 


(c) For working fluid vapor dimensionless velocity: 


ie oe eee. 0 vV*¥ = 0.0 
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BB. AUXILIARY EOUATHOnNS 

A means to derive saturated working fluid heat of vapori- 
Zation, density, and pressure is required before a numerical 
solution is attempted. Using methanol as the working fluid and 
data from Eubank [Ref. 8] for saturation conditions, the 
following functions are derived. 

1. Polynomial for the Working Fluid Heat of VaporiZzatpieg 

Tabulated values of the heat of vaporization for methanol 

are plotted on Figure 3.1 with a solid line. A Fortran program, 
polyfit, is used to obtain the coefficients of the following 


polynomial: 


He ee at eee eee (Il Ten 


where a = -9.96 x1074 


Huy AU Or) b = .730 Btu/Lbm(°R) 7 age 
c = 397.99 Btu/Lbm. The values of slag Ee using this polynoniam 
are plotted on Figure 3.1 with a dotted line. 


2. Polynomial for the Working Fluid Vapor Saturation 
Density 


Tabulated values for methanol saturated density are 
plotted on Figure 3.2 with a solid line. Again the Fortran 


program, polyfit, is used with the following modification: 





Nie T where ID ee = 500°R and.) —aec 


A polynomial in the variable € is used to describe the 


saturation density. 
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Mien polynomial, Sa(é), 1s then: 


Se (2) 16S aie? 3 eS 5 clams boc. 2A) 
Meee 202 tbm/ft-. b = 2100475 Lbm/ft~, and 
© = 2.8868 Lbm/£t> amd. GQ = —2. 7050 Lbm/ft>. The value of 


dSa/dT is also required for the numerical solution. MThis is: 


dsa _ (Zac? + Mele" + c&)/T 


aT (Tet 25.) 


ref 


eee ©Olynomial for the Working Fluid Vapor Saturation 
Pses sure 


Tabulated values for methanol saturation pressure are 
Pltected On Figure 3.3 with a solid line. Polyfit is used with 
the same modification. 


The polynomial Fa(&) 1s: 


Fa(&) = Bee + eo ele. 4: 
weer 4.5994 Lbf/in-, B= -191.2472 Lb£/in~, 
Meee G40 Tbt/in-, and @ = 439.9883 Lbf/in*. ‘The first 


derivative of Fa(&) with respect to temperature is equal 


eo: 


(3a ey a5 2bE- ae (allay (Gleiaies 2 6) 


ref 
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IVa UESCUSS TON ANDT RECOMMENDATIONS 


pee DISCUSSION 

The final version of the model for a variable conductance 
heat pipe with a vertical orientation is a system of three 
ordinary, non-linear, coupled, differential equations. These 
meee ECuations are (111.14), (172.19), and (III.22). 

The remaining challenge for the numerical model is to pro- 
pose an iteration scheme that can solve this system of three equa- 
tions. A special thank you goes to the collaboration of Dr. M.D. 
Kelleher and Dr. K.T. Yang for the following iteration technique. 

The steps of the iteration procedure are as follows: 

(a) To start the iteration process, assume a form for the 


relationship between 6 and Ws such that a C18 + C.- thas 


linear relationship is reasonable because: 


5 om When n = O, both We and 6 are values close to 
zero. 

(11) When n = L both Wo and 96 are equal to 1.0. 

(111) This step is required for the first iteration 
Omily . 


The equations for this step are then: 


= C19 r Co Crys) 
dw 

a _ dé 
“dn ad Cy dn ° CE 2 } 
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Then, using Equation (111.19) , solvemten = This establishes 


a first degree equation, in terms of 68 only which Can @emem 


do 
dn’ 
be integrated numerically to establish a first iteration 
relationship for "Giipe 

(b) With this relationship for 9 and _ defined from step 
(a) above, Equation (111.22) can be se@vedmeec. ae 

(c) Once step (b) is completed, Equation (II1.14) cangspe 
solved for dw /dn. The integral of dw /dn is then used to 
define W- The process can then be repeated. 

In summary, the first iteration goes through steps a-c 
and the remaining iterations go through steps b-c until there 
are no further changes in the form of 6, Woe and V*s ain 
theory, the outcome of this numerical iteration is a table 
which gives dimensionless temperature, 6; velocity, V*; and 


mass fraction Ofsthe working fluid vapern, Woe versus non- 


dimensional condenser length, n. 


B. RECOMMENDATIONS 


1. This numerical analysis was written to yield results via 
a computer program. The Continuous System Modeling 
Program III [Ref. 9] is a computer package that could 
be used to solve the system of equations. The CSMP 
[Ref. 9] includes a powerful integral statement which 
carries out a Fourth Order Runge-Kutta Method. 


2. This model assumes constant values for the binary mass 
diffusion coefficient and the convection heat transfer 
coefficient. Once numerical values are obtained from 
the iteration scheme, the accuracy can be improved by 
using analytical expressions instead of constants for 
the binary mass diffusion coefficient [Ref. 10] and the 
convection heat transfer coefficient Rerae ae 


3. An attempt could be made to model the variable conduc- 


tance heat pipe that is not ina vertical orientation. 
Kelleher [Ref. 3] shows that a stratification of the 
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working fluid and noncondensible gas occurs when the 
heat pipe is in the horizontal position. This position 
makes the transport processes inside the pipe three 
dimensional, which would require a finite difference 
computer model for solution. 


on 
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